function meanR2 = model_CalMeanR2(model)
meanR2 = sum(model.meanAlpha2+model.meanBeta2)/(2*pi);